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Abstract 

Hilbert Spaces of bounded one dimensional non-linear oscillators are stud- 
ied. It is shown that the eigenvalue structure of all such oscillators have the 
same general form. They are dependent only on the ground state energy of 
the system and a single functional X(H) of the Hamiltonian H whose form 
depends explicitly on H. It is also found that the Hilbert Space of the non- 
linear oscillator is unitarily inequivalent to the Hilbert Space of the simple 
harmonic oscillator, providing an explicit example of Haag's Theorem. A 
number operator for the nonlinear oscillator is constructed and the general 
form of the partition function and average energy of an non-linear oscillator 
in contact with a heat bath is determined. Connection with the WKB result 
in the semi-classical limit is made. This analysis is then applied to the specific 
case of the x 4 anharmonic oscillator. 
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I. INTRODUCTION 



In this paper we shall study the general structure of the energy eigenvalues for one 
dimensional non-linear oscillators. To be specific, we are interested in Hamiltonians which 
have the form, 



H 



eJala + ^J +V{a,a^) , (1) 



where V(a,a^) is the interaction Hamiltonian and is a functional of a and a\ the creation 
and annihilation operators for the simple harmonic oscillator (SHO), and eo is the SHO 
energy scale. We shall restrict ourselves to bounding potentials for which V(x) — > oo when 
\x\ — > oo. When V is a polynomial in a and at consisting of terms (a)) p a q the degree I of V 
is the maximum value of I = p + q for the polynomial. As is well known, when I < 2 the 
Hamiltonian is easily diagonalizable by either shifting the operator by a constant (for 1 = 1) 
or by a Bogoluibov transformation (for 1 = 2). Nonetheless, the results of these analysis 
have had far reaching applications, including coherent and squeezed quantum states [0J in 
quantum optics, and the theories of superfluidity and superconductivity 0. 

When I > 2 the oscillator is usually called anharmonic with the classical example being 
the oscillator with an x A or (a + a^) 4 interaction potential. This particular non-linear oscil- 
lator has been extensively studied since the early 1970's ( 0- 0; see |7J for a review of 
the literature), due mainly to the equivalence between it and the 4 quantum field theory in 
one-dimension. It is hoped that a detail study of this simplified system will shed some light 
on the structure of 4 theory in higher dimensions. Research on this oscillator continues 
today, mainly because it provides a natural test bed for such approximation schemes as 



the strong coupling expansion |§, modified perturbations schemes, 0, [fLOl , variational 



modified perturbation theories |IT|], lattice methods |12[], etc. More recently, Bender and 



Bettencourt have provided a deeper understanding of the system by using multiple- 



scale perturbation theory showing that the frequency of oscillation depends on the energy 
H of the state. This was interpreted by them as an operator form of mass renormalization. 
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The main purpose of this paper is not to present a new method of calculating the energy 
eigenvalues of non-linear oscillators, although we shall end up doing so. Rather, it is to study 
the general structure of both the Hilbert Space and the energy eigenvalues of non-linear 
oscillators with arbitrary binding V(a, at). The approach we shall take follows most closely 
the analysis done for the I < 2 oscillators. Namely, we shall attempt to construct, in much 
the same way, operators a and a) from a and at which diagonalizes the Hamiltonian. We 
find that unlike the SHO operators, a and at obey the commutation relation [a, at] = X(H) 
where in general A (if) is a functional of H . Its precise form depends on the specific choice of 
H and is a constant only when I < 2. The study of any non-linear oscillator thereby reduces 
to the study of operators having this commutation relation along with the determination of 
X(H) and the groundstate energy of the system. 

Because A(i?) is not a constant function in general, we find that a and at cannot be 
unitarily equivalent to a and at . Only in the special case when I < 2 does such a trans- 
formation exists. Consequently, the Hilbert space of the non-linear oscillator is generally 
unitarily inequivalent to that of the simple harmonic oscillator. This is an explicit example 
of Haag's Theorem, first proposed by Haag in 1955 |TJ[] (see also [|I5[]) for quantum field 
theories. In this paper Haag actually proved a weaker version of the theorem by showing 
that the unitary transformation between the non-interacting and interacting quantum field 
theories via the interaction picture does not exist. Later, this result was extended by Hall 
and Wightman [|16| (see also [fL7| ) who showed that based on the Wightman axioms the 



expectation values of the product of four or fewer fields of an interacting theory are unitarily 
inequivalent to those of the free theory. 

There have been other attempts at using algebraic methods to analyze non-linear os- 
cillators, of course, such as the action angle or time operator methods (see for example 
fl8|j - [p2]|). Both of these methods, however, are generalization of classical analytical tech- 
niques to quantum mechanical systems. They rely on the existence of the phase <fi and time 
T operators which are canonical to the number and Hamiltonian operators: [N, (/)] = 1, 
[H, T] = 1. Because of the positivity of the spectrum of both N and H for bound systems, 



such operators do not exist in the usual quantum mechanical system |23| (see, however, 
p4| - PH[| for the existence of such operators in extended quantum mechanical systems). In 
this sense, these methods of solution are "formal". The approach we have taken in this 
paper does not suffer from these problems. It is not a generalization of classical techniques 
but is instead a generalization of the Bogoluibov transformation and is inherently quantum 
mechanical in nature. Classical solution techniques such as the action angle are used only 
in the semi-classical limit where they are expected to be valid. 

The rest of this paper is organized as follows. In Sec II the general Hilbert Space and 
energy eigenvalue structure of non-linear oscillators are analyzed. It is found that both 
depend on a functional A (if) of the Hamiltonian. A number operator is constructed and 
the Heisenberg equations of motion are solved. Then in Sec III thermal or KMS states are 
analyzed and it is shown that both the partition function Z and and average energy (H)t 
for non-linear oscillators are similar in form to those of the SHO. In Sec IV a method of 
determining A (if) is outlined and in Sec V the connection between this method and the 
semi-classical WKB result is shown. Application of this analysis to the x 4 potential can 
then be found in Sec VI. Concluding remarks are given in Sec VII. 

II. GENERAL STRUCTURE 

Given a Hamiltonian H constructed from a, and a\ we seek solutions of the operator 
equation 

[a,H] = e X(H)a, (2) 

where a is understood to be a functional of a and . This is an eigenvalue equation with 
a being the "eigenoperator" of H and A(if) its corresponding "left eigenvalue", although 
unlike the standard eigenvalue equation A (if) is a functional of H and the ordering in eq. (Q) 
is important. Eq. (fj) does not determine a uniquely since if a satisfies eq. (Q), then so does 
g(H)a and ag(H) where g is any functional of H . A normalization for a is needed which we 
choose to be 
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H = e {a) a + e s ) , (3) 

since it diagonalizes the Hamiltonian explicitly. eoe g is the ground state energy of the 
system and is a constant. This is very similar to the way one determines the Bogoluibov 
transformation which diagonalizes the I = 2 Hamiltonian H = e^a + iei(a 2 — (a^) 2 )/2, but 
now A (if) is a functional of H . With this normalization, eq. (Q) reduces to 

[a,a f ] = \{H). (4) 

To show that a and a) creates and annihilates eigenstates of H, we make use of the 
identity, 

[a, H n ] = {(\(H) + H/e ) n - (H/e ) n }a . (5) 

obtained using eq. (£|). Then for any given functional M(H) which is expandable in a Taylor 
series, 

[ti, M(H)] = {M(X(H) + H/e ) - M(H))~a . (6) 

With this we see that if a is an eigenoperator of H with left eigenvalue A(if), then so is a n , 

[a n , H] = eS\{H) + \{\{H) + H/e ) + ■■■ + 

A(A(- ■ ■ A(A(F) + H/e ) ■■■ + H/e ) + F/e Q )|a" . (7) 

Given eq. (0), the Hilbert Space TC n i for non-linear oscillators and their energy eigenvalues 
are easily constructed in much the same way as the SHO Hilbert Space H-sho- Namely, if 
\4>)ni is an eigenstate of H, then so is a n \4>) n i as well as (a^) n |^») n /. Since the spectrum of the 
operator a) a must be non-negative it is straightforward to show (see Appendix A) that if 
A is a positive definite function then there exists a state (the groundstate) \VL) n i in TCsho 
for which a\Q) n i = and has energy eoe g |2l|. H n i is therefore spanned by the states 

\n)nl - 7^= , (P) 



where 

A n = X(e g ) ■ (x(e g ) + A(A(e g ) + e g ) 

X(e g ) + A(A(e g ) + e g ) + ■ ■ ■ + A(A(- • • A(A(e s ) + e g ) ■ ■ ■ + e g ) + e g )J . (9) 
They are eigenstates of H with eigenvalues eoe n where 
e n = e n _i + A(e n _i) , 

= e g + X(e g ) + X(X(e g ) + e g ) + ■ ■ ■ + A(A(- • ■ A(A(e 9 ) + e g ) ■ ■ ■ + e g ) + e g ) . (10) 

A thereby determines the splitting between successive energy levels. 

If X(H) is a constant, then from eq. ([II]) we see that the energy levels of the oscillator are 
equally spaced. As we shall show in the next section, this is only possible for I < 2, which is 
well known. When I > 2, X(H) is a functional of H and this equal spacing no longer occurs. 

Notice, however, that both eq. @ and the commutation relation eq. (§|) are invariant 
under unitary unitary transformations: a — > UaW. As usual, unitary transformations are 
canonical transformations which preserves the commutation relation. For the SHO, A = 1, 
while for an non-linear oscillator X(H) is a functional of H. Since a unitary transformation 
cannot change the functional form of A, a and a are unitarily inequivalent. Consequently, 
the Hilbert Spaces T~C n i and Tisno ar & unitarily inequivalent Hilbert Spaces. 

It is well known that under certain conditions any two solutions of the canonical com- 
mutation relation [a, a'] = 1 is connected by an unitary transformation. Our result does not 
contradict this. The operators a and which create and annihilate eigenstates of a general 
one dimensional oscillator do not in general obey the canonical commutation relation except 
in the special case of / < 2. Instead, [a, a)] = X(H). Indeed, it is precisely for this reason 
that a and a cannot be related to one another by an unitary transformation. 

If we consider now Hamiltonians of the form eq. (|l|) in which V(a,a^) is controlled by 
a single coupling constant e 1; we can label the Hilbert Space for each e\ as 7i ei . Then 
7i tl is unitarily inequivalent to TisHO- Moreover, there cannot be a unitary transformation 
which maps 7i ei — > Ti^ when e\ ^ e[. If there were, then using a succession of these 
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transformations we can construct a unitary transformation mapping 7i ei to TisHO and the 
two Hilbert Spaces would be unitarily equivalent. Thus for different values of e\ the Hilbert 
Spaces 7i ei are inequivalent to one another. This is a concrete example of Haag's Theorem 
(see [[HJ), first proved for quantum field theories using Lorentz invariance. 

Because the spectrum of cJa no longer consists of the non-negative integers, but instead 
depends on the energy of the state, a)a cannot in general be interpreted as the number 
operator N for non-linear oscillators. We shall construct this operator explicitly by first 
noting that H does not change the occupation number n of the states |n) n j. Consequently, 
[N,H] = with the subsidiary condition that N\Q) n i = 0. N(H) can be a functional of H 
only. From eq. (§), it is straightforward to see that it must satisfy the commutation relation 

[a,N(H)] = a. (11) 

Expanding N(H) in a power series and once again using eq. the solution to the operator 
equation eq. (PD reduces to finding the solution of the algebraic equation 

N(X(e) + e) — N(e) = 1 , (12) 

for a given A with the "boundary condition" N(e g ) = (see Appendix B). e is a real 
number in eq. fljjD and N(H) is obtained by replacing e^e — > H. (Or, equivalently, eq. fll2|) 
is the resultant equation after applying the corresponding operator equation to an eigenstate 
of H with energy e e.) 



Like the differential equation it resembles, the general solution of eq. ([12]) consists of the 



linear combination N p (H) + Nh(H) where N p (e) is the "particular" solution to eq. (|12"|) while 
Nh(e) is the solution to the "homogeneous" equation 

N h (X(e) + e) - N h (e) = . (13) 

Whether or not such a solution exists depends on the particular form of A(e), although the 
above equations can be solved for general A in the semi-classical limit as we shall see in 
Sec. V. Unlike a differential equation, however, the single boundary condition N(e g ) = 



is not sufficient to determine N(H) uniquely in general. Consider the case of the SHO. 
Then A = 1 and the particular solution of eq. ( |P2D gives N P (H) = H/e — 1/2 and is, in 
fact, the usual number operator. However, the solution of the homogeneous eq. ( p~3| ) is any 
periodic function with period 1 which vanishes at e = e g . There are an infinite number of 
such functions, such as 

N(H) = H/e -e g I + C sm(Tt(H/e - e g I)) . (14) 

for any real C. 

From eqs. and (§) the study of non-linear oscillators reduces to the determination 
of the groundstate energy e g and the functional A (if). This is non-trivial and a method 
for doing so will be given in the Sec IV. For now we shall limit ourselves to a qualitative 
description of the energy levels by looking at different possible behaviors of A(e). 

For a ground state to exist, A(e 9 ) > and we shall restrict our considerations to such A. 
Representing the eigenvalues of H/eo generically by e, if A(e) is a monotonically increasing 
function which is unbounded from above, then the energy spacings between successive energy 
levels becomes wider as n increases and e n grows rapidly with n. If, on the other hand, 
A(e) — > constant as e — > oo, then eventually the energy levels become equally spaced and we 
would once again obtain SHO type of energy levels. Notice also that if we consider eq. ( fLOf) as 
a non-linear transformation of e„ generated by A(e), then the fix point of this transformation 
A(A(e) + e) = A(e) occurs precisely when A(e) goes to a constant (see Appendix B). Finally, 
if A(e) is a monotonically decreasing function of e which decreases sufficiently rapidly, there 
will be an upper bound to the energy levels e max . 

Finally, let us consider time evolution. If H does not explicitly depend on t, time 



evolution is generated by a unitary transformation pi . 



a(t) = e itH/h a{0)e- itH/h , (15) 

which preserves the commutation relation eq. (^). Using eq. (Q), the solution to the Heisen- 
berg equation of motion is 
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a(t) = e- ieoX ^ H)t/h d(0) . (16) 

The frequency of oscillation of a(t), eo\(H)/h, now depends on the Hamiltonian H. This 
agrees with the recent result of Bender and Bettencourt |]I3|| and was interpreted by them 



as an operator form of mass renormalization. A n 's dependence in eq. (H) on the energy of 
the state would then be wavefunction renormalization. 

III. KMS STATES 

We now put the non-linear oscillator in contact with a thermal reservoir at a temperature 
T and consider the average energy and number density of the system. We shall denote these 
thermal averages by (■ • -)t which we shall take to be a KMS state [|3(J- |3^| . Namely, if 



A(t) and B(t) are two operators in the Heisenberg representation, then 

(A(t)B(t)) T = (B{t)A(t + ih/3)) T , (17) 
where 1/(3 — ksT '. Since thermal equilibrium states are stationary, we can without a loss of 



generality take t = in eq. ( |T7l) 



Applying this condition to eq. (|J), 

(X(H)) T = (a(0)a t (0)) T - (a(0)a\ifi(3)} T . (18) 

Then using the solution eq. ( |Ifi| ) of the Heisenberg equation of motion, eq. (0) and the 
commutation relation eq. (f|), 

(\{H)) T = ((H/e - e g ) (l - e'^ H % , (19) 

which reduces to the usual Bose-Einstein distribution for the SHO when A = 1. 

Unlike the case of the SHO it is not possible to determine (H)t any further using solely 
the KMS condition. We must make use of a partition function and shall restrict ourselves 
to states which can be represented by a trace over a density matrix, 

(H)t = \^n nl He~? H , (20) 



where Z = Tr^ n; e ^ H is the usual partition function. Then using the identity 

([H + e X(H)]e-^ x ^) T = {e~^ x ^) T {H) T - ^L( e -«W0) T , (21) 



we obtain 

(H) T = e e g - l _ {e _ eomH))T , (22) 
and we see once again the importance of A (if). Indeed, from eq. ( p0|) we find that 

e -/3<E e 9 

Z = 1 _ ( e - £O 0A(tf)) T • ( 23 ) 

As for the number operator, from eq. (0), 

(e- eo/3AW ) T = (iV(Jf + e X(H))e~^ x{H) ) T - (N(H) e - e °^ H) ) T . (24) 

Then using 

(N(H + e \(H))e-^ x W) T = ±Tr Hnl N(H + e,X{H))e^ H '^ x ^ , (25) 
and eq. (gDJ), we find that (N(H + e X(H))e- eol3X ^} T = (N(H)) T so that 

( e -«oA(H)) T = (l - e -^W) ) t . (26) 

This once again agrees with the SHO result for A = 1. 

IV. SOLUTION OF THE EIGENVALUE PROBLEM 

a and A (if) can be determined in the following manner. Since if is given in terms of a 
and a\ in general a = a(a, a'), which is understood in terms of a power series, 

oo 

a= Yl b rs (a?) r a s . (27) 

r,s=0 

By using the commutation relation [a, a'] = 1, we can always reduce any expansion of a to 
this form. Eq. ( p?D is well defined only if the corresponding function 

10 



oc 



f(z,z)= E b rsZ r Z% (28) 
r,s=0 

is absolutely convergent on R 2 . 

At this point we should also express A(if) as a power series in H/eo, insert this series and 
as well as eq. (p7|) in eq. (0) and obtain an infinite set of coupled equations between various 
b rs and the coefficients of the A expansion. The problem would quickly become intractable, 
however. We shall therefore first make the following drastic simplification. Instead of eq. (0) 
we shall solve the simpler equation 

: [a, H] := e \(H)a , (29) 

where :: denotes normal ordering. Correspondingly, we shall take the normalization condition 

as 

: H := e (: ~a)a : +~) . (30) 

We shall then use the solution of this equation as a guide to reconstructing the solution to 
eq. (H). Notice that corrections to the groundstate energy cannot be determined under this 
simplification and can only be determined when the full operators are reconstructed from 
the solution to eq. ([29]) . 

Denoting the solution to eq. fl29"|) by the superscript sc, we find that for 

CO 

g (-c) = bi s s c \a)) r a s , (31) 

r,s=0 

we have 



eo E bifX{H){a^ r a s = £ 6^ c) e (s - r)(o + )V + s : (o^'o'-^a, V{a, o + )] : + 

r,s=0 r,s=0 I 

r : (a t ) r ~ 1 [a t ,\/(a,a t )]a s : | . (32) 
Under this normal ordering, solving eq. (|32D is equivalent to solving the differential equation, 
A<«. (e<«>) /<».(,, ,) = {e W /M}PB s 8^«gl _ , (33) 
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where a^ sc ' = f( sc \a,a^). is obtained from H/cq by replacing everywhere a — > z and 
a* — > 2 and the normalization condition (|30|) is now 

e («0 = |y( sc )|2 + 1 / 2 . The right hand 



side of eq. ( p3|) is just the classical Poisson bracket but with the generalized coordinates 

where m is the mass of the particle. We are therefore looking for a semi-classical solution to 
eq. (H). Indeed, we shall see explicitly in Sec. V that the solution of eq. (|33|) is equivalent 
to the WKB approximation. 

Importantly, eq. (|3~3|) has the same symmetry properties as eq. (|2]). Namely, if /( sc )(z, z) 
is a solution to eq. (|33|), then so is f^ sc '(z,z)g(e^ sc ') where g{eS sc ^) is any function of e^ sc ' 
(although they do not satisfy the normalization condition eq. (|30D). Making use of this 
symmetry, we change coordinates to and 9 = —ilog(z/z)/2 from \z\ and 9. Then 
eq. ( |33"D reduces to 

iA W(eW)/W( e W, S ) = U^, (35) 



whose solution is 



/(-)( e M 0) = JeM-l/2exp jf f ^yj # (36) 

and satisfies the normalization condition explicitly. Determination of a {sc) is then reduced to 
performing the above integral, which requires inverting the equation = e^ sc '(\z\, 9) and 
solving for \z\ in terms of and 9. Next, for /( sc ) to be analytic on R 2 , /( sc ) fe^,0^ = 
f(sc) (V^tt), giving 

1 1 f 2n /(9e (sc) V 1 

A(-)(e(-)) = 2^ Jo \d\zf) # ' (37) 

which determines \( sc \ Notice that in contrast to phase angle techniques which require the 

construction of a phase operator (see, for example |[22|| ) and its concombinant difficulties, 

our analysis uses the phase only in the semi-classical limit where it is well defined. 

Reconstruction of a and X(H) from f( sc \z, z) is now straightforward, although tedious. 

can be obtained by first expanding /( sc ) in eq. (|36|) in a power series in z and z, then 
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taking d^ sc ^ = f( sc '(a, a)). Since dS sc ^ was obtained via normal ordering, there is an ordering 
ambiguity when we reconstruct a from it. Fundamentally, this arises when we replace z — > a, 
z — > in /( sc ^(a, a*) since the term zz in the expansion can be replaced by either a^a or 
aa^ = 1 + a) a. Therefore, to determine a we shall take , but we now 

I. J order 

replace zz — > a) a + A where the A's are constants. These are determined by requiring that 
the resulting expansions for a and A (if) satisfy both eq. (Q) and eq. @ (or equivalently 
eq. (^|)) term by term in the expansion. This uniquely determines not only a and A (if), but 
e g as well. 

From eq. ( |3?| ) we see that for A(if) to be independent of H, de^ /d\z\ 2 = k, where k 
is a function of 9 only. This limits I < 2. Correspondingly, if / > 2, A (if) is necessarily a 
functional of H. 



V. THE WKB APPROXIMATION 



We now make the connection between the solution of eq. (|33|) and the semi-classical 
limit. From the correspondence principle, in the large n limit e n goes over to the classical 
result. The spacings between energy levels e n — e n _i are small in comparison to e n _i and 
the levels are essentially continuous. In this limit, we can then approximate 



de 
dn 



where e(n) considered as a continuous function of n. Then from eq. ([ToD 

de 



dn 



A (sc) (e 



n + Boo ~ — 



I re r2n ( Q e ( sc )' 



de {sc) d6 



(38) 



(39) 



where we have replaced A — > in this limit. Integrating and using eqs. (0) and 



(40) 



2n Jo Jo \ d\z\ 2 j 

where is an integration constant which can be neglected in the limit n — > oo. Changing 
variables back to z and z in the integrand of eq. ([40]) and using eq. (^4|) , we find that 



13 



n + noo^— / / dxdp. (41) 
2nn J Jv eQe 

The integration is now over a disk T> eoe centered about the origin in the classical phase space. 
This is just the semi-classical Bohr- Sommerf eld quantization rule obtained from the WKB 
approximation. 

We next consider the solution of the algebraic eq. (|T2| ) in the large n limit. Since A(e) 
measures the energy splitting between energy levels, in this limit A(e) C e (see Appendix 
A) and eq. (|S|) can be approximated by the differential equation 

A (sc) (e)^^l, (42) 

where once again we have replaced A — > \( sc \ The solution to this equation is trivial and 
we once again obtain the WKB result, 

N(e) +71^^— : [ [ dxdp . (43) 
Inn J Jv toe 



Notice, however, that now the operator N(H) can now be obtained directly from eq. (^3| 
by expanding the integral in powers of e and replacing e — > H/e . 
Finally, we consider the quantum partition function 



Z = Tr Hnl e-P H = £ e"^" , (44) 

n=0 

in the large temperature limit. Making use of the Euler-Maclaurin formula, 

/•CO 1 . . 

J e~ l3eoe dn+ -e~^° + (e"^ ) . (45) 
In the large temperature limit ksT 3> eoe g we can neglect the terms ~ e _/3e °. Moreover, at 



this energy scale fc^T, e n 3> A(e n ). Then using eq. (|39|) , we convert the integral over n to 



one over e. Making use once again of eq. (|34|) , we find that 



Z « / / e^ eoe dxdp, (46) 
2tt?i J J 

where the integral is over the classical phase space. This is precisely the classical result with 
the requisite factor of the fundamental phase space volume 2irh. 

14 



VI. THE X A INTERACTION 



In this section we shall apply the above analysis to a non-trivial system: the x A anhar- 
monic oscillator, 



H = ^(ata+^ + ^a + at) 4 , (47) 



which corresponds to 



e (sc) = - + \z\ 2 + ^\z\ A cos A 9 . (48) 
2 eo 



Then 

iri 1(9) 



/<«) ( e <->,«) = fi^Tji^lBJ^L-j , ( 49) 



where 



m = L 7tt?5?? - (50) 

and £ = 16(e( sc ) — l/2)ei/e . When e\ > 0, this integral can be reduced to 

'<»> = Kriom F{alq) (51) 

where F(a\q) is the elliptical integral of the first kind and 



'v^K-tan 2 ^ 

yrn + tan 2 ^^ 



a = arccos . 5— , (52) 

\ / 1 i r" i -J n Z/j/ 



while 



v / T+7-l ,,,, 
9= ' (M) 



is its modulus. The analyticity of {e^ sc \9^ gives 

2 K(Jq) ' 



where K(^fq) is the quarter period of F(a\q). 
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Although our analysis is strictly valid only when t\ > 0, it is instructive to see what 
happens for t\ < 0. For |£| < 1, 

1 



1(9) 



i + viei 



F(a'\q') 



(55) 



where now 



a' = arctan 



tan 9 ^ 



and 



2v^ 



(56) 



i + vie 



(57) 



Then 



71 



(58) 



2 K(^) ' 
which vanishes when |£| = 1. 

When |£| > 1, 1(9) is complex and /^ sc ^ no longer satisfies the normalization condition 
eq. fl30l) . \( sc ^ is ill-defined. Consequently, the energy states are bounded, as is well known, 
by e<l/2 + e /(16|e 1 |). 

To determine a, we expand eq. ( flSD to third order in \z\ 3 , 

z) = z +i| {-3(, 2 - f)z + (z + zf} 

+ I ( e J.Y (^5 + 39 4 _ 25 2 3 _ U -3 Z 2 _ 3 4 
2 Ve / 12 4 8 8 

We then replace z — > a and z — > a' in the above and take 



1 



(59) 



4e 2 Ve o y 



(60) 



where 



F = -3(a 2 - (a f ) 2 )a + (a + a*) 3 + fra + f 2 a ] , 

G = ^a 5 + ?ata 4 - ^(at) 2 a 3 - 12( a t) 3 a 2 - |(at) 4 a + V) 5 
2 4 8 8 4 

+9ia 3 + g 2 a ] a 2 + g 3 (a r ) 2 a + g^f + g 5 a + g 6 at . 



(61) 
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The constants fx, f2, gi — g$ are present due to the ordering ambiguity. Requiring that 
eq. fl60P satisfies eq. (0) gives f\ — f2 — 3, while 



135 3 -153 27 

gi + gz = -15 , g% = — g- , 5-4 = -g , g5 = — g— , #6 = — ^ ■ ( 62 J 

The groundstate energy of the oscillator is also determined to this order, 

1 3d 21 /eA 2 
e 9=o + 7--^ - • ( 63 ) 



2 4e 8 \e J 
Using now the commutation relation eq. (HI), we obtain 



»(S)(MHi)'{?(! + 



69 /i? A 2 9 /# A 15' 



2 Veo 2/ 2 



(64) 



while 



3#i + £ 3 = y , (65) 

giving = 75/4 and g% = —135/4. This last relationship was obtained by requiring that A 
is a function of H only. To this order then, 

a = a + 



4e \ 


- 3 (a 


1 fei N 




2 W 


>"{§ 


75 , 
4 


135 


8 



ata 2 -^(at) 2 a-^(at)3_l^ a _| a t|. (66) 



4 

The energy levels can now be straightforwardly calculated from eq. (|10D . 

e n = e n _i + A(e n _i) , 



n-l 



e 



+ E K*) ■ (67) 



r=0 



Using eq. (CT), and keeping terms to order (ei/e ) 2 only, we obtain after re- arrangement, 



ei n_1 

e n = e g +n + 3—^2(r + 1) 

e o r=0 



eA 2 /51^ 2 51^ 21\ , , 
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Then 

13ei /; , x (ti\ 2 {17 , 51 , 59 21\ , . 

e n = n + - + z -(2n* + 2n + 1) - (-) (^— + ^ + y n + -) , (69) 

which is the standard second order perturbation result. Notice also that if we keep terms 
only up to ei/eo, then e n — e n -\ = A(e n ) ~ 1 + 3ne\/eo. This is precisely the result obtained 
by Bender and Bettencourt |T3l . 



It is instructive to compare eq. (plj) with the expansion of eq. (pl|) 



A^(e) = 1 + 3 (^) (e< sc > - 1/2) - ^ f ^ V ( e W - 1/2) 2 . (70) 



Notice that in both expansions the coefficients of the highest power of the energy in each 
term are the same. This is a generic feature. Quantum mechanical corrections to A ( - sc - ) only 
results in the appearance of lower powers of H/eo in each term of the expansion. Moreover, 
if we then use A^ to calculate e n , we find 

e ( n sc) =n+\ + ^(n 2 - n) - (*\* ( ^n 3 - f n 2 - ^ , (71) 



2 2e v ' \e J V4 
which also agrees with eq. (57) in the large n limit. This also is a generic feature of the 



expansion since the coefficient of the highest power of n in each term of the expansion is 
obtained from only. 

The above perturbative result is valid only for small e± and n. In the large n limit, the 
semi-classical result is valid and 

A ^" AM ^"^(vf- (72) 

To compare with the WKB result, from [|7| we know that 

WKB 34/3?r2 ( £A „4/3 f7Q\ 

e " [r(i/4)f/ 3 U • (73) 

This gives the energy splitting between levels as 

ApWKB a 3/2 / WKB, \V4 
gWJB _ g WKB „ tl « £2 !i . (74) 

n+1 3 n [r(l/4)] 2 V eo / 

Since K(l/\/2) = |T(l/4)] 2 /A/ y/n, this is precisely the form of X(e n ) for large n and we see 
explicitly the equivalence between A ( - sc - ) and the WKB approximation. 



VII. DISCUSSION 



We have shown that the study of non-linear oscillators is equivalent to the study of 
algebras satisfying eqs. (||) and (|]); the SHO being a special case of this algebra. In addition, 
the Hilbert Space eq. (||) and eigenvalues eq. (|T0|) of these algebras all have the same form. 
The number operator for non-linear oscillators was also constructed. Results of this general 
analysis were used to determine the general form of the partition function and average energy 
for an non-linear oscillator in contact with a heat bath. Analysis of non-linear oscillators 
therefore reduces to determining the function \(H) and the groundstate e g of the oscillator. 
This can be done by first making a semi-classical approximation, which requires only the 
evaluation of a single integral, and then using it as a guide to constructing a and a) in terms 
of a and a) . This analysis was applied to the x 4 interaction and both the standard second 
order perturbation result as will as the WKB result were obtained. Moreover, the recent 
results of Bender and Bettencourt were also obtained within this framework. 

Unlike the Bogoluibov transformation, the mapping between (a, a 1 ) and (a, a)) is non- 
linear and cannot be generated by a simple unitary transformation. The two Hilbert Spaces 
H ei and TLsho are unitarily inequivalent. Indeed, each value of e± determine separate Hilbert 
Spaces all of whom are inequivalent to one another. This result provides a concrete example 
of Haag's Theorem proved first for quantum field theories in higher dimensions. Based on the 
results of this theorem and the generality of our analysis, we expect a similar construction to 
hold for the 4 theory in higher dimensions. Notice, however, that this construction requires 
a natural energy scale to define \{H). For Hamiltonians of the form eq. ([!]) we have such an 
energy scale: eo- For quantum field theories, however, no such natural energy scale exists. 
An energy scale would have to be introduced, providing a natural introduction of a high (or 
low) energy cutoff in the theory. 

Appendix A 

The proof that a ground state exists for the non-linear oscillator when A(e) > for e > e g 
follows in the same manner as that for the SHO. Let H = H/e — e g = a) a. Then, if \<p ) is 
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an eigenstate of H with eigenvalue <p , 

o = (0o|^|0o) = |a|0o)| 2 >O, (75) 

and all the eigenvalues of H are non-negative. Next, consider the state = a\<f) ). Then 

#|0_ 1 > = 0o|0-i>-A(#)|0_ 1 >. (76) 

Since A is a functional of H + e g , from the above must be an eigenvalue of H also, 

which we shall label as </>_i. Then eq. ( |76| ) reduces to 0_i = O — A(0_i +e 9 ). Since A(e) > 0, 
we have 4>-i < <A> 

Similarly, the states \<f>- n = a n |0o) are also eigenstates of ^ with eigenvalues 0_ n . More- 
over, they satisfy a sequence of strict inequalities 

0_ n < 0_( n -l) < • ■ ■ < 0-1 < </>o • (77) 

Since 0_ n > for all n, this sequence must end. Namely, for some m, <p m = 0. Then 

2 

= (0_ m |i?|0_ TO ) = a|0_ m ) > or a\(j)- m ) = 0. The ground state is then identified as 
|fi) = |0„ m )and H\n) = e e g \n). 

The existence of a groundstate is only gaurenteed when A is positive definite. In the 
semi-classical limit it can be shown that this holds for wide classes of bounding potentials. 
For the polynomial interaction potential V(a, a*) — 6i(a + a*) 1 /I, 

A W ~ e 1 ' 2 - 1 ' 1 ( e A 1/l . (78) 



For the exponential interaction potential V(a,aJ) = e 1 e a2( - a+a ^ 2 



In both cases A^^ is positive definite. Notice also that A^ c -* / e^ sc ' — ► as — > oo, justifying 
the approximations made in SecV. 

Appendix B 
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The form of eq. @ as well as the algebraic equation eq. (|i~2|) suggests that we look at 
a finite difference form of the one dimensional Lie derivative. Given any function X(x), we 
define the finite difference Lie operator C x by 

C x f(x) = f(X(x)+x)-f(x), (80) 

where f(x) is any function of x. In the limit X(x) — > for all x, it is straightforward to see 
that eq. (|80| ) reduces to the usual Lie derivative. Moreover, for any two functions f(x) and 
g(x) and constants a, b, 

C x (af(x) + bg{x)) = aC x f{x) + bC x g(x) . (81) 

£\ is therefore a linear operator on the space of functions on R 1 . It is not, however, a 
derivation since it does not satisfy the product rule, 

C x {f(x)g(x)} = f(x)C x g(x) + g(x)C x f{x) + C x f{x)C x g{x) . (82) 

Finally, for any two given functions X(x) and £(x), the commutator of two finite difference 
Lie operators 

[C x , Cs\f{x) = f(x + X(x) + e(x + X(x))) - f(x + £(x) + X(x + £(x))) , (83) 

vanishes if and only if C x (,(x) = C^X(x) . 

Using C x , the particular solution of eq. ([12]) becomes the solution of the operator equa- 
tion, 

C x N p = 1 . (84) 

Of more interest is the homogeneous solution to eq. (|T2|). 

C x N h = , (85) 

which lies in the kernel of C x , ker C x , for a given A. Notice that when A is a real constant, 
ker£^ is the space of all periodic functions with period A. When A is a function of x, ker/^ 
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will contain generalizations of periodic functions to those whose frequencies are x dependent. 
This agrees quite well with the observation that in the semi-classical limit, l/X^ reduces 
to the WKB result. Of particular interest is when A G ker L\. 

£\\ = . (86) 

From eq. (|lOl) we see that for this A the energy levels are equally spaced and is determined 
solely by A(e g ). The constant function, and thus the SHO, clearly satisfies eq. (|86|) . Whether 
there exists other non-trivial solutions to eq. ( [36"1) for physically realizable non-linear oscil- 
lators is still an open question. 
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